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ABSTRACT 



Aims. We have investigated the properties of the central star and dust in the bipolar nebula IRAS 19312+1950, which is an unusual 
object showing characteristics of a supergiant, a young stellar object, and an asymptotic giant branch (AGB) star. 
Methods. We obtained //-band polarimetric data of IRAS 19312+1950 using the near-infrared camera (CIAO) on the 8 m Subaru 
telescope. In order to investigate the physical properties of the central star and the nebula, we performed dust radiative transfer 
modeling and compared the model results with the observed spectral energy distributions (SEDs), the radial profiles of the total 
intensity image, and the fraction of linear polarization map. 

Results. The total intensity image shows a nearly spherical core with ~3" radius, an S-shaped arm extending ~ 10" in the northwest 
to southeast direction, and an extended lobe towards the southwest. The polarization map shows a centro-symmetric vector alignment 
in almost the entire nebula and low polarizations along the S-shaped arm. These results suggest that the nebula is accompanied by a 
central star, and the S-shaped arm has a physically ring-like structure. From our radiative transfer modeling, we estimated the stellar 
temperature, the bolometric luminosity, and the current mass-loss rate to be 2 800 K, 7 000 Lq, and 5.3 x 10~ 6 MQyr~', respectively. 
Conclusions. Taking into account previous observational results, such as the detection of SiO maser emissions and silicate absorption 
feature in the 10 yum spectrum, our dust radiative transfer analysis based on our near-infrared imaging polarimetry suggests that (1) the 
central star of IRAS 19312+1950 is likely to be an oxygen-rich, dust-enshrouded AGB star and (2) most of the circumstellar material 
originates from other sources (e.g. ambient dark clouds) rather than as a result of mass loss from the central star. 

Key words, near-infrared - imaging polarimetry - radiation transfer modeling - individual (IRAS 19312+1950) - (post-)AGB stars 
- circumstellar matter 



1. Introduction 

IRAS 19312+1950 is a peculiar bipolar nebula with a horn- 
like appearance extending ~30" in the 2MASS JHK-band im- 
ages. With IRAS colors [12] - [25] = 0.50 and [25] - [60] = 
0.71, it falls into region VIII on the IRAS color-color diagram 
( Ivan der Veen & Habingll988h . A relatively large fraction of ob- 
jects in this region is known to belong to different types such 
as Miras, OH/IR stars, or carbon stars. A weak 10 yt/m absorp- 
tion was detected in the I RAS Low-Spectr al Resolution spec- 
tra of IRAS 19312+1950 dVolk et alj Il99lb. Maser emissions 
of sili con monoxide (SiO), water (H2O) (Nakashima & Deguchi 
2000), and hydroxyl (OH) (Lewis, B. M., private communica- 
tion) were also detected. These observational results led the au- 
thors to conclude that IRAS 19312+1950 is a candidate for an 
oxygen-rich asymptotic giant branch (AGB) or post-AGB star 
(Nakashima & Deguchi 2000). However, some unusual charac- 
teristics have been found in subsequent observations, and de- 
bates on the t rue na ture of IRAS 19312+1950 have arisen. 
Deguchi efal] d2004l) found C-, N- and O-bearing molecules 
(e.g. H 13 CN, CH3OH, SO and S0 2 ) towards this object. C- and 
N-bearing molecules such as H 13 CN and CH3OH are often de- 
tected in carbon-rich circumstellar environments or in molec- 
ular cloud cores in a dark cloud. While C-bearing molecules 
and SO take ~10 5 yrs and ~10 5 yrs, respectively, to be formed 
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in an environment like a dark cloud, the time scale (~10 yrs) 
of expanding gas with an expansion velocity of 10 km s _1 to 
cross th e envelope of IRAS 1 9312+1950 (~5.8xl0 17 cm at 
2.5 kpc, Naka shima et al.ll2004l) is remarkably shorter than the 
above time scales. Furthermore, a stellar velocity of ~36 km s _1 
was found from the velocity components of the SiO and H2O 
masers. Besides the broad component in the thermal emission 
of some molecules, the presence of a narrow component in 
the thermal emission of CO, HCN, CS, SO and HCO + sug- 
gests the presence of an additional, kinematically cold compo- 
nent, which has been interpreted to originate from a dust cloud 
rather than the matter ej ected by mass loss fro m the central star 
dNakashima & Deguchill2000t iDeguchi et al.|[2004b . In addition, 
iDeguchi et al.l (2004) estimated the envelope mass and the mass- 
loss rate to be ~25 Mq and ~2.6 x 10~ 4 M0yr _1 , respectively, 
by means of the large-velocity gradient (LVG) model with their 
observational results of several molecular lines in the radio fre- 
quency domain. For the above reasons, to date, it seems plausi- 
ble that IRAS 19312+1950 is a progenitor with an initial mass 
of M* > 4 M(?) and is embedded in an ambient cloud by chance 
(see also Nakashima & Deguchi 20051). 

Previous studies of this object have mainly focused on identi- 
fications of gas-phase molecules and investigations of their kine- 
matics. In order to provide a better understanding of this object's 
nature, model analyses to explain observed spectral energy dis- 
tributions (SEDs) and images as well as to investigate the stel- 
lar parameters are important. Although dust is expected to con- 
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tribut e only ~ 1 % of the mass in the envelope (e.g Knapp & Kerrl 
1 19741) . it plays a more important role in the transfer of radiation 
from the optical to far-infrared (FIR) than gas. Furthermore, light 
scattered by dust in the envelope often dominates in the SED at 
optical and near-infrared (NIR) wavelengths, and the maximum 
dust grain size in envelopes around AGB stars an d post-AGB 
stars is expected to be on the order of 0.1 /mi (e.g. lJural [l996). 
Polarization analysis in the NIR is a powerful technique for 
probing such a dusty environment. In our experiments, therefore, 
we performed imaging polarimetry of the IRAS 19312+1950 
nebula in the NIR and radiative transfer modeling to derive the 
stellar parameters and dust properties in the envelope by com- 
paring these with the observed SEDs, the total intensity image, 
and the polarization properties. The results of our observations 
and numerical simulations are presented in sect.|2]and[3] respec- 
tively. In sect. [4] we discuss the properties of dust in the envelope 
and the central star of IRAS 19312+1950. 



2. Observations and Results 

2.1. H-band imaging polarimetry 

We obtained polarimetric images of IRAS 19312+1950 using 
the NIR camera (CIAO) on the 8 m Subaru telescope on 2 
August 2002. Because the object is faint in the /-band (mj = 
11.3 mag) and the nebulosities are easier to detect in the H band 
than the /f-band, we only tried the //-band in this experiment. 
The central wavelength and the band width of the //-band fil- 
ter are A c = 1 .65 //m and AA = 0.30 /mi, respectively. We used 
the medium resolution camera with a 0.0217 arcsec pixT 1 pixel 
scale (22" x 22" field of view (FOV)). Since the target is too 
faint at the wavelength for adaptive optics (AO) wavefront sens- 
ing (m«=18 mag), the AO was not used. The natural seeing was 
~Cf.'6 at //-band. We f ollowed the observing sequence described 
in our previous paper (Mura kawa et al.ll2005l) . In order to mea- 
sure linear polarization, we obtained four sets of images with 
orientations of the half-wave plate (HWP) of 0°, 45°, 22.5°, and 
67.5°. The exposure time for each frame acquisition was 3 sec. 
Five frames were obtained for each orientation of the HWP and 
eight dither offsets of 10" separation were performed. The total 
integration time was 480 sec. The nebula of the target covers the 
entire FOV of the camera, and co mplex interstellar dif fuse emis- 
sion and polarization are expected dAxson & Ellisl 19761) because 
the target is located near the galactic plane (galactic coordinate 
of 55°. 37, +0°.185). Thus, we also obtained sky frames for cali- 
bration of the sky background level and interstellar polarization. 
The position offsets are 20" east and west of the target. The same 
observing procedures as those of the target were applied. 

We reduced the observed da ta in a similar manner as de - 
scribed in our previous papers ( Mur akawa et all l200l l2005h . 
We first subtracted the dark frame and applied flat field- 
ing. We calibrated the sky background level and the interstel- 
lar polarization by subtracting the sky frames. Then we ob- 
tained the Stokes IQU images, converting with formula / = 
(h.o + hi.s + /45.0 + hi .5) 12, Q = (/o.o - /45.0) In, and U = 
(hi.s - hi. 5) fa- Where 77 is the total polarization efficiency 
of the HWP and the wire-grid polarizer, the value of 0.875 
dMurakawa et al.l 120041) is used in this experiment. The polar- 
ized intensity (PI), the fraction of linear polarization (P), and 
the polarization position angle (ff) are also derived by PI = 
Ve 2 + U 2 , P = PI /I, and 6 = l/2arctan(t//g), respectively. 
We also obtained error images of all Stokes parameters, the 
fraction of linear polarization, and the polarization position an- 
gle. The signal-to-noise ratios of the Stokes / image are 7-10 



in a region within 1 " from the central star and 30-60 in a re- 
gion ~3" from the central star. The estimated background sur- 
face brightness and the polarization in the neighboring sky are 
7 x 10~ 15 Wrrr 2 /fln _1 arcsec -2 (13.0 mag arcse c" 2 , cf. 13.4 mag 
arcsec -2 for the averaged sky background level, Tokunaga 2000) 
and 1-2 %, respectively. The error of linear polarization is 2-5 % 
per pixel. 



2.2. Polarization map 

Figure Q] (a) and (b) present a total intensity map (Stokes /) 
with polarization vectors and a fraction of linear polarization 
map, respectively. The field of view (FOV) is 15"xl5". Figure 
|2] shows what is seen in our images. The total intensity map 
shows a nearly spherical nebulosity within 3" of the inten- 
sity peak, two ring-like structures extending ~10" to the north- 
west and the southeast (PA of -37° with respect to the in- 
tensity peak), and a lobe extending southwest. The ring-like 
structures look rather point-symmetric. These nebulous fea- 
tures are seen as a horn-like structure in t he 2MASS JHK- 
band images (Nakashi ma & Deguchil l2000h and a NIR im- 
age dDeguchi et al. .1 12004 . The polarization vectors are centro- 
symmetrically aligned in the entire nebula, where a sufficient 
signal is detected, and the fraction of linear polarization attains 
30-60 % in nebulosities at NW, SE and SW. These features sug- 
gest that the nebulosity is seen with singly scattered light in the 
//-band, which is emitted from the central star, and the nebu- 
lae are associated with the central star. In the fraction of linear 
polarization map, we see low polarization (P < 20 %) in the 
SE and NW arms. A possible interpretation is scattering angle 
dependence of polarization properties. The fraction of linear po- 
larization has the maximum value near the normal scattering an- 
gle, and low values in forwa rdly or backwardly scattered light 
(see Fig. 2 in the paper by iFischer et alj[l994f) . If we assume 
that the S-shaped arm has a physically ring-like dust structure 
and the plane of the ring is tilted with respect to the plane of 
the sky, the scattered light from the ring is expected to have 
low polarizations. Such a polarizat ion distribution is remar kably 
different from IRAS 17441-2441 dOppenheimer et alj|2005l) . in 
which polarization is enhanced at the rim of the bipolar lobes. 
In this object, the bipolar lobes are hollow and light scattered at 
the rim in the plane of the sky results in a large fraction of po- 
larization. The different characteristics in the polarization map 
of IRAS 19312+1950 are due to a different density structure of 
its nebula. We also identified two spots with high polarization 
within 1" of the central star. The spots are aligned in the equato- 
rial direction; i.e., perpendicular to the polar direction. A similar 
signature is seen in an optically thi n nebula with an equatoria lly 
enhanced density distribution (e.g. Muraka wa et alJlin prep.l) . It 
may be possible that dust condenses more in the equatorial di- 
rection than in the polar axis. However, it is slightly suspicious 
because of the natural seeing-limited angular resolution (~076 at 
//-band) in our image. 

We compare our NIR image with the re sults of previous 
CO observations (Naka shima & D eguchi 2005). Narrow compo- 
nents of 35-36 and 37-38 km s -1 13 CO J = 1-0 emission lines 
extending ~ 20" roughly correspond to the NW and SE arms. 
Because the velocity distribution does not suggest any rotating 
motion, the narrow components could represent a bipolar out- 
flow rather than a disk or a spherically expanding shell. The NW 
and SE arms are blue- and red-shifted components, respectively, 
while the NW arm is fainter than the SE one. It can not be ex- 
plained with the inclination effect of the nebula. Broad compo- 
nents in 12 CO and 13 CO only appear around the central star with 
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Fig. 1. Results of //-band imaging polarimetry of IRAS 19312+1950. The field of view is 15" x 15". a) The total intensity (Stokes /) 
image. The polarization vectors are also plotted every 14 pixels (~073). b) The fraction of linear polarization map. The poor image 
quality 2" from the bottom is due to the lack of data in some dither positions. 
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Table 1. Model parameters 



Fig. 2. This cartoon explains what is detected in our images. 



~ 5" radii. It could reflect the inner core, as seen in our NIR 
image. We will discuss more details of these narrow and broad 
components in sect. 14. 21 



3. Radiative transfer modeling 

3.1. Model assumptions 

We performed radiative transfer calculations to estimate the stel- 
lar parameters and to investigate the dust properties in the en- 
velope. We used a Monte Carlo code, which simulates ther- 
mal emiss ion and light scatte ring by dust grains in a model 
geometry dOhnaka et alj 12.006b . This code computes the SEDs 
and a dust temperature distribution, and generates Stokes IQUV 



parameters 


adopted values 


comments 


stellar model 


T e s 


2 800K 


by SED fit 




7 000 L Q 


by SED fit 


d 


2.5 kpc 


adopted 


R* 


2.5xl0 13 cm 


calculated 


a* 


0.66 mas 


calculated 


model geometry 




1R, 


by SED fit 


R u 


50 R* 


by SED fit 


Rout 


23 000 R+ 


adopted 


T 


20 (V), 1.5 (H) 


by SED fit 


A 


2 


adopted 


Pi 


0.8 


by SED fit 


grain model 


^min 


0.005 yum 


adopted 




0.3 yum 


by polarization fit 


n(o) 


a- 35 


adopted 



images. Our code was applied for interpretation of the mid- 
infrared spectro-inte rferometric data of the silicate carbon star 
IRAS 08002-3803 dOhnaka et alJ ?20061 and the near-infrared 
polarimetri c data of the bipolar prq to-planetary nebula (PPN) 
Frosty Leo (Murak awa et alJlin prep.l) . 

Although attempts to estimat e the stellar parameters have 
been made, they are still uncertain. Nakashi ma et alj {2004) have 
estimated a distance ranging between 2.5 and 3.9 kpc by apply- 
ing the bolometric luminosity of 8 000 Lq, which is typical for 
AGB stars. The lower limit is obtained if the estimated interstel- 
lar extinction of Ay = 16.3 mag from their CO integrated inten- 
sity (Visr = 27-33 km s component) is assumed to be from 
a foreground cloud. If it is fr om a background cloud , the upper 
limit is obtained. Following Nakashi ma et al.l {2004), we apply 
the distance d of 2.5 kpc in our modeling and discuss the effect 
of the distance on the luminosity in sect. l4.1I The stellar effective 
temperature is also not well known, but is expected to be less 
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than 4 000 K, taking into account the possible spectral type of 
a late K or an M giant (Wood, P. R. 2002, private communica- 
tion). The bolometric luminosity and the stellar temperature are 
estimated by SED fitting. 

The detection of seve ral oxygen-rich molecular lin es (e.g. 
masers of SiO and H2O, INakashima & Deguchil |2000|) and a 
weak absorption feature at 10 pm dVolk et aTll 1991b indicate an 
oxygen-rich environment in the envelope. Therefore, we assume 
bare silic ate grains. W e use the optical constants of astronomical 
silicate (Draine 1985). To simplify the grain model, we applied 
a spherical shape with an MRN-like power law size distribution; 
i.e., a m j n < a < a max and n (a) = a -3 5 , where a and n (a) are 
the grain radius and the size distr ibution function, respectively 
(Mathis, Rumpl & Nordsieckl 19771) . In our Monte Carlo experi- 
ments, the size distribution averaged values of the opacities and 
the scattering matrix components are used, which are given by 

< X >— I X(a)n(a)da I n(a)da, 

where X is the extinction cross section C ext , the absorption cross 
section C a b s , or the scattering amplitude matrix elements S 1 (9) 
or S 2 (#)■ The averaged scattering cross section C sca is given by 
< C sca >=<Cext>-<C'abs>- We set the minimum grain size a m { n to 
be 0.005 pm since there is no crucial clue to determine it. The 
maximum grain size a max is constrained by our modeling. 

The IRAS 19312+1950 nebula already shows a distinct and 
complex bipolar appearance in the ~ 10" scale FOV. However, 
in our modeling, we focus on estimating fundamental properties 
such as the stellar temperature, the bolometric luminosity, the 
dust envelope mass, and the present-day mass-loss rate instead 
of reproducing the ring-like structure. Because the inner core 
(r < 3") has a nearly spherical structure, although an indication 
of asymmetry is expected, we apply one-dimensional spherically 
symmetric geometries, and such an approximation would satisfy 
our science goal. The outer radius R out is set to be a fixed value 
of 23 000 R*, which corresponds to the extension of the nebula 
of ~15" detected in the 2MASS JHK-band images, assuming 
the distance of 2.5 kpc. Since two origins of dust (i.e., the stel- 
lar system and ambient clouds) have been suspected before, we 
assume a model geometry with a radial density profile given by 
p (r) = p\ (r/7?i n ) ^' + pi (r/Ri n y^ 2 , where is the inner radius 
of the dust shell and the densities of two components become 
equal at a transition radius R tr . The values of pi and pi are deter- 
mined from an input value of the radial optical depth of the shell 
(r). 

3.2. Model results 

We found that exponent values of fii =2 and /fe = 0.8 pro- 
duce the best fit of the mid-infrared SED, in particular the 
self-absorption silicate feature around 10 pm. With respect to 
the grain properties, we obtained the maximum grain size of 
flmax = 0.3 (+0.05) pm, which satisfies the fraction of linear po- 
larization of our observational result. The extinction cross sec- 
tion Cext, the dust albedo to, and the scattering asymmetry g- 
parameter at 1.65 pm are 8.47 x 10~ 7 cm 2 , 0.389, and 0.187, 
respectively. We searched the following parameter ranges: T e g 
of 2 500 to 3 000 K with a 100 K step, U of 5 000 to 10 000 L Q 
with a 1 000 L© step, R m of 3, 5, 7, and 10 R*, R tl of 10, 30, 
50, and 100 and ry of 7, 10, 20, and 50. These combinations 
yield 2 304 models, and we have examined their SEDs. The best 
fit model has r eff = 2 800 (+100) K, U = 7 000 (+1 000) L©, 
t v = 20 (1.5 at 1.65 pm), R in = 7 and/? tr = 50 The stel- 
lar radius R+ and its apparent size at 2.5 kpc are 2.5 xlO 13 cm and 
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Fig. 3. A comparison of the model SED with the observational 
results. The solid, dotted, dashed, and dashed-dotted curves de- 
note the model results of the total flux, scattered light by dust in 
the envelope, the attenuated stellar flux, and the thermal emis- 
sion from the envelope, respectively. The filled circles, filled 
squares, open crosses, and filled triangles are from the photo- 
metric results using SIRIUS on the UH88 telescope and MSX, 
the IRAS low resolution spectrum (LRS), and the IRAS pho- 
tometry, respectively. The estimated interstellar extinction of 
Ay = 16.3 mag is corrected in the SIRIUS data. The esti- 
mated background contributions of 30 % and 62 % of the 60 and 
100 pm flux are subtracted from the IRAS photometry. 



0.66 mas, respectively. The model parameters are summarized in 
TableQ] 

The SED of the best model was compared with the ob- 
served SEDs in Fig. [3] The data are from from JHK-band pho- 
tometry obtained using the SIRIUS camera on the University 
of Hawai'i 88-inch telescope (UH88), which is dereddened 
with Ay = 16.3 mag t o correct for the foreground extinction 
(INakashima et al.ll2004 . the Mid-course Space eXplorer (MSX) 
point source catalogue, the InfraRed Astronomical Satellite 
(IRAS) low-resolution spectrum (LRS), and the IRAS photom- 
etry. Because IRAS 19312+1950 is located near the Galactic 
plane, a large frac tion of flux in FIR is expected from back- 
ground emission (Izumiuraet al. 1999). We obtain a rough es- 
timation of the background contributions of 30 % at 60 pm and 
62 % at 100 pm from the IRAS Sky Survey Atlas (ISSA). Since 
the pixel scale of the images is 115, the target is unresolved and 
seems to be located in a region with complex emission features; 
the real values could exceed the above values. These background 
contributions are subtracted from the results of 60 and 100 pm 
IRAS photometry, plotted in Fig. [3] Scattered light from the dust 
shell and the stellar flux dominates at wavelengths shorter than 
5 pm, while the thermal emission from the dust shell does at 
longer wavelengths. The dust temperatures at the transition ra- 
dius (50 /?*) and at the inner boundary of the shell (7 7?*) are 
290 K and 960 K, respectively. The inner boundary tempera- 
ture suggests a steady formation of fresh dust in the outflow as 
oxygen-rich dust starts to condense at temperatures typically in 
the range 1 000-1 500 K. The flux from the inner shell mainly 
contributes at wavelengths shorter than 10 pm and the flux from 
the outer shell at longer wavelengths. Our model SED is fair 
agreement with the IRAS L RS 10 pm spec trum, which shows 
the silicate self-absorption dVolk et alj 1 199lh . The FIR flux of 
our model results are sufficiently lower than the background 
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(a) 



model result 




1 2 

separation (arcsec) 




1 2 

separation (arcsec) 



Fig. 4. The radial profiles of the total intensity (a) and linear 
polarization (b), respectively. The solid curve and the dashed 
curve denote the observational result and the model result, re- 
spectively. 



subtracted results of the IRAS photometry at 60 and 100 /jm. 
The possible reasons are that our model geometry (i.e., one- 
dimensional spherical model) is too simple and the estimation 
of background contribution is too low, as already described. 

The applied distance of 2.5 kpc in the above analysis is based 
on an assumption that all CO emissions of Vi sr — 27 - 33 km s" 1 
component are from foreground clouds (Nakashima et al. 2004). 
It is hard to determine the distance to the object in general. 
However, it may be possible to be constrained by taking the 
contribution of the interstellar extinction into account. If all CO 
emissions are assumed to be from bac kground, the distance t o 
the object is estimated to be 3.9 kpc (Nakashi ma et alj 12004). 
This results in all foreground extinctions are due to the interstel- 
lar extinction o f Ay = 6.2 mag applying the stan dard interstellar 
extinction (e. g. lAllen|[l97l IDeguchi et al.ll 1 998b and the optical 
depth of the envelope of t> = 31. We performed some mod- 
eling with different contributions of the foreground extinction 
between 6.2 and 16.3 mag. We found that a model with Ty > 23 
(i.e., a foreground extinction Ay < 14.9 mag) fits the observed 
SED, and a large Ay is preferred more than a small Ay. Thus, 
it is likely that the object is near the back side edge or behind 
ambient clouds. Then we take into account the contribution of 
the interstellar extinction with a fixed foreground cloud compo- 
nent. We found good solutions in SED fit when the interstellar 



extinction Ay is between 2.2 and 6.7 mag. From this result, we 
estimate a possible distance to be between 1 .4 and 4.2 kpc. 

We also made Stokes IQU images with another Monte Carlo 
experiment. Since the thermal emission is negligible in the H- 
band, as seen in the SED (Fig. [3), we simulate with monochro- 
matic light at a wavelength of 1 .65 fim to obtain Stokes images 
at a particular wavelength much more efficiently than the full ra- 
diation transfer mode, which is applied to calculate SEDs. In this 
fixed wavelength mode, a photon packet is always scattered at an 
interacting point, and the weight of the photon packet is reduced 
by a factor equal to the dust albedo co. The output Stokes im- 
ages were convolved with a Gaussian function with a 076 diam- 
eter, which corresponds to the observed PSF size in the //-band. 
Figure [4] (a) shows a comparison of radial intensity (Stokes /) 
distribution of our model to the azimuthally averaged observa- 
tional result. Our observational result can be actually fit with 
/3 = 1 .2 if a single power law is applied. However, it provides 
a poorer fit of the 10 /urn silicate self-absorption feature and too 
low flux in the FIR and we could not find any adequate solu- 
tions to fit the observed SED and the radial intensity distribu- 
tion simultaneously with a single power law model. Although we 
see somewhat steep and shallower slopes in our two component 
model within 079 from the central star and further, respectively, 
the two component model can reproduce the observed data rea- 
sonably well. Figure|4](b) shows the radial profile of the fraction 
of linear polarization. Since there are regions with small polar- 
izations along the S-shaped arm, we omitted the data towards the 
central star within r < 173 because of contamination by unpolar- 
ized light from the central star PSF. On the other hand, the linear 
polarization gently increases for r > 173. Although our observa- 
tional result of the radial profile has a somewhat complex shape 
because of the S-shaped arm in the bipolar nebula, the overall 
shape of our model result is in good agreement with our obser- 
vational data. Changing the maximum grain size does not affect 
the shape of the radial profile much, but it does change the ab- 
solute fraction of linear polarization dramatically. For a smaller 
«max < 0.25 /mi, a larger polarization > 45 % is attained, and for 
a larger a max > 0.37 /im, a lower polarization < 35 % at r — 475 
is attained. 

4. Discussion 

The puzzling thing about IRAS 19312+1950 is that this ob- 
ject shows characteristi cs of several object classes. The re- 
ported massive envelope (Deguchi et al. 2004) and an extremely 
high mass-loss rate are reminiscent of a supergiant. However, 
some C- and N-bearing molecul es such as NH3, SO, H2CS, 
and CH3OH (IDeguchi et al.ll2004l) in the envelope are ones of- 
ten detected in dark clouds. The detection of SiO (and H2O) 
masers and a stronger OH maser satellite line at 1612 MHz 
than main lines at 1665 and 1667 MHz suggest an AGB star 
or a post-AGB star. Furthermore, an FIR emission distribution 
(> 1') extends sufficiently more than the nebulosity seen in 
the NIR (~ 30")- As already suspected, these characteristics 
are probably caused by contributions from t he matter originated 
outside the object; i.e amb ient clouds (Deg uchi et alj 12004; 
Nakashima & Deguchi 2005). Thus, we will discuss the prop- 
erties of the central star and the dust in the envelope separately. 

4. 1 . Stellar properties 

From our radiative transfer modeling, the stellar temperature 
and the distance are estimated to be T e ff = 2 800 K and d = 
1.4 - 4.2 kpc, respectively. The uncertainty of the distance 
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causes the estimated luminosity (1 10CW 2 Lq) to range between 
2 200L© an d 19000Lq- With these properties, the object class 
of IRAS 19312+1950 could be classified in the Hertzsprung- 
Russell diagram. Considerin g YSOs, low mass YSOs are much 
fainter (L* < 10 Lq, e.g. iHartigan et~ai1ll995t IWhitnev et ail 
120041) . and luminous YSOs such as Herbig Ae/Be stars are 
much hotter (T e ff > 4000 K) (see also Fig. 1 of theoreti- 
cal pre-main sequen ce tracks in the HR diagram presented by 
|Palla&Stahlerll999l) . Thus, IRAS 193 12+ 1950 is unlikely to be 
a YSO. Considering evolved stars, it is also unlikely to be a su- 
pergiant because the estimated luminosity of I RAS 19312+1950 
is within the AG B limit (~ 5 x 10 4 Lq, iPaczvnskil Il97lt 
IWood et aT1ll983l) and is sufficiently lower t han several know n 
supergiants; e.g. 7 x 10 5 L q for IRC +10 4 20 dJones et al.ll 19931) . 
2 x 1 s Lq for NML Cy g dBlockd|200ll), 5.5 x 10 5 Lq for VY 
CMa (Sopka et al. 1985). Therefore, the possibility remains that 
IRAS 19312+1950 is an AGB or a post- AGB star. 

From our modeling, we estimate the dust envelope mass 
to be 0.14 Mq if the distance to the object of 2.5 kpc is ap- 
plied. Comparing th is with a gas envelope mass of ~25 Mq 
(Deguchi et al. 2004), the gas-to-dust mass ratio is found to 
be ~ 180 : 1, which is very typical for the neutral interstel- 
lar medium in the galactic plane and the wind s of AGB star s 
of about solar metallicity dKnapp & Kerjll974t lKnapplll985l) . 
We also estimate a mass-lo ss rate. While the v alue of 2.6 x 
10 -4 Mq yr _1 obtained by iDeguchi et alj d2004l) is derived at 
r ~ 3" from the central star, our modeling allows us to de- 
rive the current value. Applying the above gas-to-dust mass 
ratio and the expansion velocity of 25 km s , which is ap- 
proximately half t he width of the CO J = 1-0 emission ve- 
locity distribution (Nakashi ma et al J 12004). we obtain the cur- 
rent mass-loss rate of 5.3 x 10~ 6 Moyr -1 . This value is in 
good agreement with the value of 6.3 X 10~ 6 MQyr -1 calculated 
with an empirical formula of mass-loss rate logM = -5.65 + 
1 .05 l og(L/10000Ln)-6 3 log(r eff /3 500 K) for M-type AGB 
stars dvan Loon et alj|20"05l) . From the above properties, the cen- 
tral star of IRAS 19312+1950 shows characteristics of dust- 
enshrouded AGB stars. 



4.2. Nebula and duality of the dust chemistry 

Dust shells formed in the AGB phase often have a nearly 
spherical shape, as seen in TTCyg (Olofs son etaT] [2000) 
and in large field of view (> 10") images of IRC +10216 
(Mauron & Huggins 2000). Bipolarity, which is seen in a large 
fraction of PNs dCorrad i & Schwarz 1995), is thought to be 
form ed by the fast-wind during the PPN phases (e.g. review 
bv lBarick&Frankll2002l) . However, the evidence for the pres- 
ence of high velocity jets and a rotating disk have been reported 
in some AGB stars, which are not thought to undergo the in- 
tensi ve mass loss that forms an aspherical dust shell; e.g., V 
Hya dKahane et al.ll 19961: ISahai et al.ll2003b iHirano et al.ll2004l) . 
X Her dKahane & Jural [19961: iNakashimal f2005h . and RV Boo 
(Bergman et al.l |2000|) . These facts are hard to explain with a 
scenario of stellar evolution of a single star, and it is thought 
that binary companion(s) must play an imp ortant role in shap- 
ing aspherical shells (e.g. ISokerlll997l 120021) . With respect to 
IRAS 19312+1950, the S-shaped arm is detected in NIR im- 
ages, and such a morphology might be a result of interactions 
with binary companions. However, while the high velocity broad 
components are extended and seen as a bipolar outflow, and nar- 
row components concentrate close to the central star in case of V 
Hya, the kinematic properties in IRAS 19312+1950 is opposed 
sense. A possible interpretation is given as follows: the object is 



probably embedded in an ambient cloud. The fast wind from the 
central star, which is detected in the broad components in 12 CO 
and 13 CO, rapidly slows down due to the thick ambient cloud 
matter. As a result, the narrow components of 35-36 and 37- 
38 km s -1 in 13 CO J = 1-0 extending ~ 15" towards the NW 
and SE directions are seen. 

Duality of dust chemistry, i.e., detection of C- and O-bearing 
molecules, in the IRAS 19312+1950 nebula is also an inter- 
esting issue. Carbon stars with oxygen- rich molecules in the 
envelope, such as silicat e carbon stars dLittle-Mareninl F l986: 
Willems & de Jondfl986l) and the Red Rectangle dWaters et all 



1998), are known or are expected to have binary companions. 



The most widely accepted hypothesis of the presence of oxygen- 
rich dust in envelopes around these carbon stars suggests that 
silicate was ejected by mass loss from the primary star and was 
trapped in a circumbinary disk or a companion disk when the 
prima ry carbon star was an M giant ( Morris ll 987b iLlovd-EvansI 
1990). Evidence for the presence of a disk or a reservoir is 
obtained with detection of a narrow velocity componen t (v < 
5 km s -1 ) in CO emission lines (e.g. Kah ane et al.lfl9 98). With 
respec t to IRAS 19312+1950, the detection of SiO maser emis- 
sions (Nakashima & Deguchi 2000) presents an important con- 
straint on the nature of the central star. SiO maser emissions 
have been detected toward evolved stars with only a few excep- 
tional star forming regi ons; e.g., O rion- KL dBaud et alJll980h 
and W5 1 IRS2 dFuente et al.ll 19891) and are pr oduced abov e the 
surface (5-10 AU or a few #*) of evolved stars dElitzurll980l) . To 
date, no detection h as been reported in a silicate car bon star (e.g. 
i Nakada et al.ll 19871) or in a ny other carbon star (e.g. lLepine et al.l 
ll978l;ISchoieretalJl2006l) . Thus, IRAS 19312+1950 is proba- 
bly an oxygen-rich star at least at this moment and C-bearing 
molecules do not originate in the mass loss from the central star. 
The origin of C- and N-bearing molecules, which are detected 
tow ar ds IRAS 19312+1950 could be chemical reactions in the 
wind dLindqvist et al.|[l988l: I Willacy & Millarlll997l;[Duari et all 
1999) or ambient cloud matter, which kinetically merges with 
the matter ejected by the mass loss from the central star. 



5. Conclusion 

We present H-band imaging polarimetry of the peculiar bipolar 
nebula IRAS 193 12+ 1950 using CIAO on the Subaru telescope. 
The total intensity image clearly shows a point-symmetric S- 
shaped arm extending towards the northwest and the southeast. 
The polarization map revealed a centro-symmetric polarization 
vector pattern in the entire bipolar lobes with polarizations of 
30-60 %. These results indicate that the nebulae are associated 
with the central star and are seen in light scattered by dust in 
the nebula. We find low polarizations (P < 20 %) on the east 
side of the SE and NW arms, suggesting that the S-shaped arm 
has a physically ring-like structure instead of a hollow or cavity 
structure in the lobes. 

We also investigated physical properties of the central star 
and the nebula by means of dust radiative transfer calculations. 
The estimated stellar temperature and bolometric luminosity are 
2 800 K and 7 000 Lq, respectively. The distance uncertainty is 
found to be between 1.4 and 4.2 kpc. Comparing the gas enve- 
lope mass of 25 Mq to the dust envelope mass of 0.14 Mq, the 
gas-to-dust mass ratio obtained is 180 : 1. With this value, the 
current mass-loss rate is estimated to be 5.3 x 10~ 6 M0yr~' . This 
is in good agreement with a value of 6.3 x 10~ 6 MQyr~' calcu- 
lated with an empirical mass-loss formula for M-type AGB stars. 
From above stellar properties and detections of SiO, H2O, and 
OH masers, the central star of IRAS 19312+1950 shows char- 
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acteristics of an oxygen-rich, dust-enshrouded AGB star rather 
than a YSO or a supergiant. 

O ur results support the idea p roposed by dDeguchi et alj 

2004; Na kashima & Deguchill2005li that the object is embedded 
in ambient clouds, because the large envelope mass of 25 Mq is 
impossible to explain with only mass loss from the possible AGB 
central star. Although the detected C- and N-bearing molecules 
can form by chemical reactions in the oxygen-rich envelope, it 
could also originate from ambient cloud matter, which kineti- 
cally merges with the outflow of the central star. 
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